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We present evidence for the absence of a gap in a class of S" = 1 antiferromagnetic 
exchange models. The spin exchange is long-ranged of the type — (— 1)'^-'/|« — j|" 
where 1 < a < 3. We have shown previously that without the alternating factor the 
model for q; = 2 (5 = 1 Haldane-Shastry model) has a gap, exponentially decaying 
correlation functions and exponentially small susceptibility at very low temperatures. 
In the case of the alternated interaction the stabilizing next nearest neighbor ferro- 
magnetic interaction changes qualitatively the behavior of the system. We have 
studied the groundstate and first excited state using a modified Lanczos algorithm 
for system sizes up to 16 sites. Also, we performed exact diagonalization for systems 
up to 8 sites and obtained the thermodynamics. The correlation functions decay with 
distance like a power law. These models define a new class of integer spin chains that 
do not show a Haldane gap. The results may be relevant to describe impurity spins 
coupled by a RKKY-interaction through a half-filled conduction-electron band. 
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Introduction. It was proposed long ago that integer and half-odd-integer spin chains 
behave qualitatively different [1]. The reason lies on a topological term in the action that 
prevails in the latter case. On general grounds half-odd-integer spin chains are gapless and 
integer spin chains show a gap [2]. This leads to important differences in the correlation 
functions. In the first case these show power law behavior and in the second case exponential 
behavior. There are however several models that do not follow this rule at special points in 
the space of interactions. In the case oi S = 1, adding a quadratic term in the interaction 
and requiring the models to be integrable, a SU{2) invariant model [3] and a SU{?>) invariant 
model (both solvable by the Bethe ansatz) [4] are gapless. Also, adding frustrating next- 
nearest-neighbor (nnn) interactions to the Heisenberg model, it has been found that for 
S = 1/2 there is a critical value acr ^ 0.2411 (where a = J2/J1, is the ratio of the nnn 
interaction to the nearest-neighbor (nn) interaction) such that for a < acr the spectrum is 
gapless (as for a = 0) while for a > acr a gap appears [5]. This has been interpreted as a 
fiuid-dimer transition. In particular, it has been shown that if the spin is half-odd-integer 
and the groundstate is translationally invariant {k = 0) the spectrum is gapless [6]. The 
dimer phase is consistent with the Lieb-Schultz-Mattis theorem [7] which states that if the 
spectrum is not gapless the groundstate should be degenerate (for half-integer spin). Also, 
recently it has been argued that translationally invariant spin chains in an applied field can 
be gapful without breaking translation symmetry, when the magnetization per spin, m, is 
such that 5 — m is an integer. It was then proposed that a Haldane gap phase can be found 
for half-integer spin [8]. It has also been shown recently that the S = 1/2-dimer chain, the 
Majumdar-Ghosh chain and the S = 1-Haldane chain are in the same phase [9]. 

Explicit tests of the validity of Haldane's proposal have concentrated on models of short- 
range interactions. Recently, we have extended this analysis to the case of long-range inter- 
actions considering the 5 = 1 Haldane-Shastry model [10,11]. Even though we might expect 
that the correlation functions should not decay faster than the interaction we found a finite 
gap and the corresponding exponential decay of the correlation functions. This may be the 
result of the frustrating nature of the interactions. The infiuence of frustration in the 5 = 1 



Heisenberg chain has also been studied recently [12]. 

The standard Haldane-Shastry model [13,14] is a periodic version of 1/r^ exchange. The 
S = 1/2 case has attracted considerable attention [15,16]. The groundstate energy and the 
correlation functions have been obtained [15,13,14] together with the thermodynamics [17]. 
The groundstate wavefunction is a spin singlet of the Jastrow-Gutzwiller form. The excita- 
tions are spin-1/2 spinous [17] that form a gas of a semionic nature [17,18]. The asymptotic 
correlations decay algebraically with exponent rj = 1 without logarithmic corrections, in 
contrast to the Heisenberg case. This indicates the absence of spin exchange between the 
spinous rendering the models solvable in greater detail than in the short-range Heisenberg 
counterpart, solvable by the traditional Bethe ansatz method. The zero-T susceptibility is 
finite [17] (and numerically the same as for the Heisenberg model) consistently with a singlet 
groundstate and a gapless spectrum. 

The S = 1 case is not integrable (for both models). In ref. 10 we studied the groundstate 
properties and the gap to the first excited state using a modified Lanczos method for small 
systems of size up to 16 spins. We obtained that the groundstate energy per spin is —1.267894 
and the value of the gap is 0.55439 (recall that for the Heisenberg model the gap has been 
estimated to be 0.41050). The general trend of the groundstate correlation functions is that 
they decay faster than those for the Heisenberg model (in the sense that the numerical values 
are smaller) both with distance for fixed N and as a function of the size of the system. The 
same happens for the S = 1/2 case where the spectrum is gapless. A linear fit of \ogC]\[/2 as 
a function of A^ yields a correlation length of the order of .^ = 3.1. This is an underestimate 
due to finite size effects. Also, we performed complete diagonalization of systems up to 8 
sites to obtain the temperature dependence of the susceptibility and correlation functions 
[11]. The susceptibility is exponentially small at low-T and the correlation functions also 
decay faster than those for the Heisenberg model as a function of temperature. 

Model. In this work we consider a class of models of long-range interactions of the type 
1/r" but with alternating signs [16]. A simple and convenient way to implement periodic 
boundary conditions is to use the chord distance (that is the distance between the points 



when the chain is wrapped into a closed circle). We consider then the Hamiltonian [19] 

where d{i — j) is the chord distance d{i — j) = \ sin(z — j)(t)\/(t) with ^ = tt/A^, where A^ is 
the number of lattice sites in the chain. The leading term (nn) is antiferromagnetic but the 
next nearest neighbor term (nnn) is ferromagnetic which tends to stabilize the dominant 
term. We take 1 < a < 3. 

This class of long-range spin exchange models may describe impurity spins in a metallic 
host interacting via a RKKY-interaction ^ cos[2kpr]/r°' if the conduction electron band is 
half-filled. The value of the exponent a in the case of non-interacting electrons \s a = d 
(where d is the dimensionality of the lattice). The case of interacting electrons has also been 
considered. In the case of the magnetic screening cloud around a single Kondo impurity 
in a Luttinger liquid {d = 1) a = Qc < I [20] and in the case of the coupling between two 
impurity spins {S = 1/2) coupled to a Id Hubbard chain a = 2 [21], for example. 

To study Model (1) we use the modified Lanczos algorithm [22] to obtain the groundstate 
properties and the gap to the first excited state and we use complete diagonalization of the 
Hamiltonian matrix to obtain the thermodynamics. In the modified Lanczos algorithm the 
size of the vectors can be considerably reduced using the symmetries of the problem. The 

-> 2 

Hamiltonian commutes with the total spin operators (St) and S^, with the translation 
operator T, the spin fiip operator R and the refiection operator L{i^N + l — i,i = 
1,...,N). The ground state has total S^ = and one of the (degenerate) first excited 
states also has S^ = 0. One can then immediately reduce the states under consideration to 
this subspace only [10]. A similar procedure can be used for the finite temperature behavior. 
We calculate the susceptibility 

X = ^E<S^.S,> (2) 



and the correlation functions 



™"iv^ 5(5 + 1) ■ ^^ 



Both diagonalizations give the exact results for the several finite-size systems. The results 
for the infinite system can be estimated using standard extrapolation methods [23] like the 
VBS method [24] or the EST method [25]. In the first method we want to estimate the limit 
of a finite sequence P„ (n = 1, ..., A^). Defining 

p(m+l) _ p{m) _| (A\ 

where Qi"^^ = 0, P^°^ = Pn we obtain an estimate of the sequence iterating. If ct^ = 
this is the Aitken-Shanks transformation which is adequate for exponential behavior. To 
generate the Pade-Shanks transformation we select a^ = 1. A power law behavior is well 
fitted choosing the Hamer and Barber's transformation a^ = ~[1 ~ (~l)'"]/2- We get an 
estimate of the asymptotic value of the sequence P„ [23,24] selecting a^ appropriately. 

In the EST algorithm we look for the limit of a sequence of the type T{h) = T + oi/i" + 
02/1^'^ + ..., where h^ = ^/N is a sequence for the several system sizes, A^. The value of the 
w}^ iteration for the sequence is obtained from 

rp(N) ^ rp{N+l) J-m-l J-m-l ,p\ 

-'m -'m-l "T / (AT+i) (AT) X \y) 

I hf, \ I 1 _ -'m-l -'m-1 \ _ 1 

UiV+mJ [^ y(iV+l) (iV+1) I i 

\ m.— l m — 2/ 

where a; is a free parameter which is adjusted such that the estimate of the error 



7^(1) _ 7^(0) I /7\ 



is a minimum (where Np is the number of data points, T_i = and Tq = T(hN))- It 
has been shown [25] that this algorithm has several advantages over the VBS algorithm in 
particular for smaller sequences, it converges faster and it is less sensitive to rounding errors. 
In this paper we will use both methods to estimate the thermodynamic limit. 

Results. In Table I we show the values for the groundstate energy and gap as a function 
of A^ (for A^ = 4 to 16 and A^ even) for the values of a = 1.01, 2, 3 [26]. In Fig. 1 we show the 
gap as a function of 1/A^ for the standard frustrating case compared to the behavior of the 



gap in the alternating case. The results suggest that these two cases are in different classes 
since the slopes for the two cases are distinct. In particular, it suggests that the alternated 
cases may extrapolate to a zero value of the gap as A^ ^ oo while the standard frustrating 
cases suggest a finite value in agreement with the results previously obtained for a = 2 [10]. 
In Table II we show the VBS estimated values for the gap using the several extrapolations. 
For small and intermediate values of a the power law fit (Hamer and Barber) consistently 
gives very small gaps suggesting that the spectrum is gapless. As a grows further (a > 2.7) 
we see that the gap increases and becomes apparently finite (we will return to this point 
later). The table shows consistently smaller gaps for the power law fit as compared to the 
other two methods. If the spectrum is indeed gapless this is consistent since, for instance 
for the Aitken-Shanks transformation, one is trying to fit a power law with an exponential. 
In Tables III-V we present the sequence of P^^ using the Hamer-Barber algorithm for 
a = 2,2.5,3. The difference between the extrapolated value and the values for the previous 
iteration gives a measure for the error involved. Due to the small number of data points 
this difference is actually an order of magnitude larger than the extrapolated value itself 
(except for a = 3). Another possible criterium is to take the difference between the values 
for m = 2, but is also of the same order. The error in this procedure is therefore large. 

To further clarify the nature of the spectrum and to have a better control on the errors 
involved we consider now the BST algorithm eq. (6). In Fig. 2 we present the values of the 
gap as a function of the free parameter u for the several values of a considered above. We 
also plot the error e defined in eq. (7). In general, for each value of a there are several points 
in (j-space where either the gap or the error go to zero or become very small. For each value 
of a there is a point (actually a narrow region) where both are very small. First we take the 
value of Lo where the gap is very small and estimate its error calculating e at that point. In 
Table VI we give the error e at the values of co chosen as above for the several values of a. 
We can also follow the standard method and select the several local minima for e and take 
the gap obtained at these points. The results are shown in Table VII. They strongly suggest 
that the spectrum is gapless (at least for a < 2.7). The BST method is consistent with the 
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VBS method for ct = 3 in the sense that it suggests a finite gap (see however ahead). 

Even though the results strongly suggest that the spectrum is gapless, it might appear 
that Table II could be consistent with a finite gap (but small compared to the frustrated 
case [10]). However, if the gap would be finite the Hamer-Barber algorithm should yield a 
finite gap as we obtained previously for the frustrated case [10] where all three algorithms 
correctly picked up the leading (finite) term of the sequence (see Table 4 of ref. 10). As a 
further check we have used the EST algorithm to find the extrapolated gap for the standard 
Haldane-Shastry model (a = 2, frustrated case). The results are shown in Fig. 3. As a 
function of co the error decreases for a; ^ 2 and the gap is saturated to the value 0.55405 
(for u = 2) very close to the (finite) result previously obtained [10]. When there is a true 
gap the extrapolation yields a finite value with a magnitude that is very close for the four 
algorithms (actually one might take the discrepancy of extrapolations of Table II as a sign 
that the gap is zero - no finite leading term in the sequence). 

In Table VIII we present the groundstate correlation functions for ct = 2. In Fig. 4a we 
show the groundstate correlation functions for a = 2 and A^ = 16 for both cases showing 
that the decay is considerably slower in the alternated case. In Fig. 4b we show the 
correlation functions C{m) for m = 1 — 4 as a function of A^ for a = 2. The behavior 
of |C(1)| is similar for the positive and alternated interactions (it is slightly smaller in 
the latter case). However, for m > 1 the correlation functions Caitifn) are considerably 
larger and reach their extrapolated values for much smaller system sizes. We estimate the 
correlation function exponent taking C{N/2) ^ (— 1)^/^1/A^'' [8] since for the system sizes 
considered it yields better results than a plot of C{m) as a function of m for fixed A^ [8]. We 
estimate r\ = 0.16,0.24,0.49, for a = 1.01,2,3, respectively (taking a fit using the system 
sizes up to A^ = 16 and excluding the A^ = 4 point). Besides the alternating signal, the 
correlation functions are modulated by an oscillatory function that decreases in amplitude 
as A^ grows. Note that \Cait{2m + 1)| > \Cait{2m)\ for m > 1 due to the oscillatory nature 
of the interaction (in Fig. 4b this is explicitly shown for |C(3)| > C(2)). 

In Fig. 5 we show the susceptibility as a function of temperature for a = 2. For A^ 



even the groundstate is a singlet but for A^ odd is degenerate (triplet) (as for the S = 
1/2 Heisenberg model). Therefore x alternates between zero and very large at small T. 
Comparing with results obtained for other models (like the Heisenberg and the Haldane- 
Shastry models for integer and half-integer spins) this alternancy suggests a gapless spectrum 
(in the case of a gap both the even and odd system sizes give a vanishing susceptibility at 
zero temperature). This would imply a finite value for the susceptibility (the extrapolation 
error is large and we do not estimate the zero-T susceptibility). As a ^ 3 we find the same 
type of behavior indicating that there is no true gap and that the extrapolated results are 
a consequence of the finiteness of the systems studied. The reason is that the decay of the 
interaction is faster (and therefore the transition to gapless behavior is slower) and the finite 
sizes considered are not enough to correctly extrapolate to zero. 

In Fig. 6 we show the correlation functions as a function of T for the positive case and 
the alternating case for a = 2. Consistently with the groundstate results the correlation 
functions also decay much slower with temperature in the alternated case. In Fig. 7 we 
show C{N/2) iov N = 8 for the several values of a. 

In summary, we have identified a new class of integer spin chains that are gapless. The 
interaction is long-ranged and non-frustrating. Indeed, the stabilizing ferromagnetic next- 
nearest-neighbor interaction changes qualitatively the behavior of the system with respect 
to an antiferromagnetic frustrating nnn interaction [10,11]. We calculated the gap (for 
5 = 1) using a modified Lanczos method finding a vanishing value in the extrapolated limit. 
We also found that the correlation functions have a considerably larger range both with 
distance and as a function of temperature. The results obtained may be useful to study 
5 = 1 impurity spins embedded in a half-filled conduction electron matrix and coupled via 
a RKKY-interaction. 
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TABLE I- Groundstate energy per spin and gap as a function of A^ for ct = 1.01, 2, 3. 

OL = 1.01 a = 1.01 a = 2 a = 2 a = 3 a = 3 

N -En/N Gap -En/N Gap -E^/N Gap 

4 2.0595832 1.1118877 2.1589760 1.2337006 2.2976817 1.3702968 
6 2.3595499 0.8899970 2.0523213 0.8480498 1.9179450 0.8541843 
8 2.5811593 0.7363615 2.0005515 0.6540099 1.7672689 0.6402738 
10 2.7580452 0.6359315 1.9707343 0.5353366 1.6920636 0.5232326 
12 2.9056771 0.5616502 1.9516741 0.4545769 1.6490979 0.4491504 
14 3.0325603 0.5042410 1.9385998 0.3957728 1.6222435 0.3979216 
16 3.1439108 0.4584054 1.9291615 0.3509020 1.6043433 0.3603571 
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TABLE II- Extrapolated values of the gap using the VBS method. 

ara = -[1 - (-l)™]/2 «™ = 1 «™ = 



a - 


= 1.01 


0.0658 


0.2004 


0.1254 


a ■■ 


= 2.0 


-0.0142 


0.1491 


0.0742 


a ■■ 


= 2.3 


0.0041 


0.1571 


0.0829 


a ■■ 


= 2.5 


0.0065 


0.1677 


0.0925 


a ■■ 


= 2.7 


0.0862 


0.1850 


0.1268 


a ■■ 


= 3.0 


0.1701 


0.2202 


0.1845 
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TABLE III- Extrapolation iterations for the gap for a = 2 using the Hamer-Barber 
algorithm. 

m 1 2 3 

1.2337006 

0.8480498 0.4575103 

0.6540099 0.3484716 0.1817683 

0.5353366 0.2825508 0.1477233 -0.0141768 

0.4545769 0.2382768 0.1244493 

0.3957728 0.2064011 

0.3509020 
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TABLE IV- Extrapolation iterations for the gap for a = 2.5 using the Hamer-Barber 
algorithm. 

m 1 2 3 

1.3002061 

0.8478680 0.4707043 

0.6421958 0.3607076 0.1988730 

0.5233554 0.2952211 0.1640630 0.0064759 

0.4452184 0.2515429 0.1406424 

0.3895432 0.2202066 

0.3476438 
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TABLE V- Extrapolation iterations for the gap for ct = 3 using the Hamer-Barber 
algorithm. 

m 1 2 3 

1.3702968 

0.8541843 0.4888593 

0.6402738 0.3818190 0.243082 

0.5232326 0.3213966 0.216709 0.1700793 

0.4491504 0.2830859 0.202198 

0.3979216 0.2570883 

0.3603571 
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TABLE VI- Values of the error e defined in eq. (7), for the several values of ct, calculated 
at the values of u where the gap goes through zero (smaller than l.OE — 06) and where e is 
also small. 

a u error (e) 

1.01 0.846448 0.00068 
2.0 0.962115 0.00017 
2.3 0.890129 0.00066 
2.5 0.768180 0.00158 
2.7 0.590684 0.00005 
3.0 0.534526 0.01318 
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TABLE VII- Minima of the error e eq. (7) for the several values of a and the corre- 
sponding values of the gap. The minima of e (smaller than l.Oi? — 05) that are close to the 
set of values of uj shown in Table VI are highlighted in bold. The absolute value of these 
results for the gap are the error estimates within the method. 

cji gap Lo-i gap (Ja gap 

a = 1.01 0.810 -0.0074257 0.864 0.0079885 1.077 0.0254297 

a = 2 0.886 -0.0067272 0.955 -0.0007323 
a = 2.3 0.733 -0.0239250 0.928 0.0034268 
a = 2.5 0.622 -0.0291923 0.883 0.0109872 
a = 2.7 0.500 -0.0163724 0.592 0.0002422 1.213 0.0713413 

tt = 3 0.914 0.1033700 1.445 0.1458258 
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TABLE VIII- Correlation functions C{m) for a = 2 for the set of values A^ = 4 to 16 
with N even. 

N 4 6 8 10 12 14 16 

m 

1 -0.75000 -0.70502 -0.68749 -0.67851 -0.67318 -0.66970 -0.66728 

2 0.50000 0.48340 0.47798 0.47560 0.47736 0.47363 0.47317 

3 -0.55677 -0.51869 -0.50316 -0.49500 -0.49008 -0.48683 

4 0.45641 0.44601 0.44094 0.43810 0.43634 

5 -0.47988 -0.46342 -0.45458 -0.44917 

6 0.43260 0.42565 0.42156 

7 -0.44604 -0.43647 

8 0.41737 
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Figure Captions 



Fig. 1- Gap as a function of 1/N for the standard frustrating interaction (positive) and 
Model (1) (alternated) for the values of o; = 1.01, 2, 3. 

Fig. 2- Gap and estimate of the error as a function of uj for Model (1) for the values of 
a = 1.01,2,2.3,2.5,2.7,3. 

Fig. 3- Gap and estimate of the error as a function of u for the standard frustrating 
interaction (Haldane-Shastry model). 

Fig. 4- a) Correlation functions C{m) as a function of m for A^ = 16 and a = 2 for 
the Haldane-Shastry model (HS) and for Model (1). b) Correlation functions C{m) for 
m = 1, 2, 3, 4 for q; = 2 as a function of N. 

Fig. 5- Susceptibility as a function of temperature for ct = 2 for the values of A^ = 

4,5,6,7,8. 

Fig. 6- Correlation functions C{m) as a function of temperature for the Haldane-Shastry 
model and for Model (1) for N = 8. 

Fig. 7- Correlation function C{N/2) as a function of temperature for A^ = 8 and a = 
1.01,2,3. 
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